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Abstract 

We compute the onium-onium scattering amplitude at fixed impact parameter in the 
framework of the perturbative QCD dipole model. Relying on conformal properties of 
the dipole cascade and of the elementary dipole-dipole scattering amplitude, we obtain 
an exact result for this onium-onium scattering amplitude, which is proven to be identical 
to the BFKL result, and which exhibits the frame invariance of the calculation. The 
asymptotic expression for this amplitude and for the dipole distribution in an onium at 
fixed impact parameter agree with previous numerical simulations. We show how it is 
possible to describe onium-e^ deep inelastic scattering in the dipole model, relying on 
/cr-factorization properties. The elementary scattering amplitudes involved in the various 
processes are computed using eikonal techniques. 
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1 Introduction 



The HERA experiments |I|, ^ have focused attention on the Bahtsky-Fadin- Kuraev-Lipatov 
(BFKL) Pomeron 0, which should be relevant for describing the small-Xf,j behaviour of the 
proton structure function. More generally, this perturbative QCD hard Pomeron describes the 
behaviour of hadronic scattering amplitudes at very high energy s and fixed momentum transfer 
t ~ — {m being an hadronic mass scale). In this Leading Logarithmic Approximation, one 
takes into account the exchange of a bound state of two reggeized gluons in t-channel. This 
resummation, based on perturbative Regge physics, predicts an increase of the amplitude 

Abfkl oc (1.1) 

where 

ap = i + ^^41n2> 1. (1.2) 
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Using the optical theorem, it follows that this behaviour violates the Froissart bound at very 
high s 

otot < c In^ s. (1.3) 

This violation is directly related to the unitarization problem of QCD, which is one of the 
main problems to be solved in the theory of strong interaction. Various approaches have been 
recently proposed in order to restore unitarity. In the multiregge approach, it has been shown 
that the Generalized Leading Logarithmic Approximation ^, ^, where one takes into account 
the exchange of any fixed number of reggeized gluons in t-channel, is equivalent to the non- 
compact Heisenberg XXX spin chain |] in the multicolor limit of QCD. The solution of this 
integrable model is still an open problem [1^, |TI|, |T2|, |TB|]. In the model recently developped 



by Mueller et al [T^, |T5|, [T6|, and separately by Nikolaev et al [|1^, |T9| , in order to control the 
perturbative approach, one deals with onia, which are heavy quark-antiquark bound states, so 
that their transverse size naturally provides an infra-red cut-off. The relevant degrees of freedom 
at high energy are then made of color dipoles. These color dipoles produce a classical cascade 
in the multicolor limit, which reveals a Pomeron type dynamics. This approach, combined with 
/cT-factorization pO|, |2l|, has been successfully applied to deep inelastic — p scattering for 
describing HERA data for F2 ^ . In the more general case of onium-onium scattering, 

the BFKL approximation corresponds to the exchange of one pair of gluons between two excited 
dipoles, each one being extracted from one of the two onia. The unitarization problem can also 
be studied in this dipole model. In real physics, unitarity and analyticity of the S matrix and 
the finite range of strong interaction lead to the Froissart bound ( |1.3D for total cross-section. 
Unitarity implies in particular at fixed impact parameter that the probability of any event 
cannot exceed 1, that is 

\Sib)\<l. (1.4) 

One way of enforcing unitarity in the framework of perturbative QCD would then be to take 
into account two and more left moving dipoles scattering off an equal number of right moving 
dipoles. As long as the relative rapidity Y of the two onia is not too high, it is possible to deal 
with unitarization effects without facing saturations effects, that is one can consider each dipole 



in the wave function of the onia as still dilute [0, ^ . However, because of non-trivial dynamics 
in transverse space, the usual Glauber multiple scattering series, which corresponds to summing 
up multiple Pomeron exchanges, diverges factorially. A possible way of escaping this problem is 
to sum over to the number of exchanged Pomerons, and only after this summation is performed. 
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averaging over configurations of dipoles in tlie onia wave functions. Such an approach, except 
for a toy model where the transverse dynamics is absent |T^, has not been yet carried out 
analytically. Numerically, such an approach is possible, and Monte Carlo simulations indeed 
show that QCD unitarizes in the dipole framework |2^, |28[| . 

Our aim is to study how this can be carried out analytically. In this paper, we will focus 
on the onium-onium cross-section at fixed impact parameter. It is organized as follows. In 
section p], we calculate the onium-onium scattering amplitude at fixed impact parameter. The 
method developped here is rather general, and is based on expansions over conformal three 
points correlation functions. Using this representation, we calculate the distribution of dipoles 
inside an onium, at fixed impact parameter, and obtain its precise dependence with respect 
to the transverse size of the excited dipole. Then, the cross-section at fixed impact parameter 
is obtained using the elementary dipole-dipole cross-section, which is calculated by eikonal 
methods and expanded on a conformal basis. Our result shows the exact equivalence between 
the dipole and the BFKL approaches for inclusive processes. In section |3|, we study some 
physical applications of the dipole picture, namely — Onium deep inelastic scattering. We 
rely on /c^-factorization and apply eikonal techniques for computing the elementary dipole-gluon 
cross-section. 



2 Onium-onium cross-section at fixed impact parameter 

In this section, we calculate the onium-onium scattering amplitude in the BFKL approximation. 
This calculation has been first carried out in Ref . ^5], |16[ . It used an expansion on a conformal 
three points correlation function basis. Asymptotic expressions of this basis were used in this 
calculation. The calculation which is presently developped gives the correct dependence with 
respect to the transverse sizes of onia. This dependence is highly non trivial and could not be 
obtained without a careful treatment of the conformal basis. 

We thus consider onium-onium scattering in the leading logarithmic approximation. The 
infinite-momentum wave function of an onium, in the large Nc limit, was calculated pertur- 



batively in Ref. [jT4[. This calculation is based on an approximation of eikonal type, due to 
the ordering of longitudinal momenta. This allows one to compute the probability of emitting 
a gluon from a quark, and then from a quark- ant iquark pair, that is from a color dipole. In 
the multicolor limit, the cascade of emitted soft gluons decouples, leading to a semi-classical 
cascade of color dipoles, in term of probability, since interference terms cancels in this limit. 
The onium-onium elastic scattering amplitude at fixed impact parameter A{Y, h) can then be 
expressed using a parton type formulation. It involves the number of dipoles in each onium 
and the elementary cross-section of two such dipoles. For a relative rapidity Y and impact 
parameter 6, AiY^h) can be expressed as 



A{Y,V} = -i J (fxij (fx2 J dzi J dz2^{xi,Zi)^{x2,Z2)F{xi,X2,Y,b). (2.5) 

^{^iLiy^i) is the square of the heavy quark-antiquark part of the onium wavefunction, x^ being 
the transverse size of the quark-antiquark pair and Zi the longitudinal momentum fraction of 
the antiquark. The momentum and P2 of the two onia are supposed to be large, with 
= = 0. Y is related to F by F = F -|- \nz1Z2, due to the fact that the perturbative dipole 
cascade originates from the quark-antiquark pair. The distribution z) of this pair cannot 
be computed perturbatively, and goes far beyond the purpose of the present approach. The 
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scattering amplitude F is then evaluated in term of the perturbative dipole cascade. Following 
Ref. |]T3|, |15], we define n(xoi, x, b, Y) such that 

Nix,b,Y) = J d\J^dzi<l>{x,zi)n{x,x!,Y,b) (2.6) 

is the number density of dipoles of transverse size x', at a tranverse distance b from the center 
of the quark-antiquark pair, where the momentum fraction of the softest of the two gluons (or 
quark or antiquark) which compose the dipole is larger or equal to e~^. Y is the relative rapidity 
with respect to the heavy quark given by y = y + ln2;i. In the leading logarithm approximation 
(noted F^^^) where the scattering is due to the exchange of a single pair of gluons between the 
two dipoles extracted from the left and right moving onia, F^^'^ reads 

F^'\x,,x„Y,b) = -l [ f-;, f-% d%, d% d'ig - 5\b, -b,- U, +g-b) 

X n(xi,x'i,Fi,6i)?2(x2,X2,F2,&2) crr>r>(^'i)^2>^i -^2)- (2-7) 

The rapidities Yi and Y2 are such that Y = Yi + Y2. Eq. ( |2.7| ) involves the elementary dipole- 
dipole cross-section at fixed impact parameter ano, which has been evaluated in [jl5|, and which 
is calculated in appendix [A.2| using eikonal techniques. For two dipoles of transverse sizes x^^ 
and Xj, whose centers are located at 6'^ and 1/2, one obtains 

f \h' h' \ ^1+^2 \\h' h' 1 1 ^ 

I |&'i-&^ + ^||&'i-&'2-^|J 

The forward scattering amplitude can then be evaluated by integration over impact parameter, 
namely 

A{Y) = f d\A{Y,b). (2.9) 



A{Y) is normalized so that the optical theorem reads 

a(r) = 2ImA (2.10) 

In order to get the expression for n{x,x', Y,b), one rehes on the global conformal invariance of 
the dipole emission kernel, related to the absence of scale. We thus decompose this distribution 
on the basis of conformaly invariant three points holomorphic and antiholomorphic correlation 
functions |2^, 3^. Introducing complex coordinates in the two-dimensional transverse space 

p={px,Py) (2.11) 

P = Px + iPy and p* = p^ - ipy, (2.12) 
the complete set of eigenfunctions -E"'^ of the dipole emission kernel is 



with 



\P10P20J VP10P20/ 



1 — n 
h = — h 

h = ——^iv (2.14) 
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being the corresponding conformal weights, with n integer and u real. This set constitutes a 
unitary irreducible representation of SL(2,C) 

We get rid of the longitudinal degrees of freedom by using a Mellin transform with respect 
to Y, namely 

n{x, x', Y,b)= J ^e'^^n^i^^ x\ h). (2.15) 
Expanding the dipole distribution on the conformal basis, one writes 



nUx,x,b)= ^ 8j—-^^lu +- n{,,„}, 



X E-''^ \^b+^-ni,b-^-wj E"'^* i^=-ni,-=-wj. (2.16) 

In this expression, the transverse integration is done with respect to the coordinate w of the 
center of mass of the quark- ant iquark pair. The longitudinal dynamics gives rise to the term 
n{u,n}uj, which was computed in Ref. and has the following expression 

2 

where 



CIril + 1 \ / Iril + 1 \ / Inl -I- 1 \ 

Ll__ + ,iy\_l^lLd_ \ =^(i)_Re^f^_ + ,z.j. (2.18) 

Going back to the longitudinal space, the distribution of dipoles takes the form 

du d'^w / 2 f^oiNc 



n{x,x,Y,b) = ^^^16j-^-^\^u+—jexpi^——x{n,u)Y 



X E-'-^ ^b+^-ni,b-^-wj E-'^* ^=-w,-=-wj. (2. 19) 

This distribution can be more easily computed by using a Fourier transform with respect to 
the impact parameter, that is by fixing the t channel exchanged momentum. We thus define 

n(x,x',r,6) = J .^e~'^-^n{x,x!,Y,q), (2.20) 



and introduce, following Ref. |2^, the corresponding mixed representation of E'^''^, namely 

^-(p) = g / f^e'l-RE-- {r + =,R-=)- (2.21) 

The normalisation term bn,u is given by 

, ^ 3 1 4.. r(-2^ + (l + |n|)/2) r{tu+\n\/2) 

-iiy+\n\/2 T{iiy + {l + \n\)/2) T{-iu + \n\/2)' ^' ' 

In this Fourier representation the dipole distribution reads 

^'.I'.B = J^/ % I ^,e-'l-tEr (x)Er(x')^e.p {^xMy) , (2.23) 



nix, X 
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2.1 The large rapidity dipole density at fixed impact parameter 

At very large Y, the term corresponding to n = dominates the exponential term in the 
expansion (|2.23|) because of the expression of i/) (see Eq. ( |2.18| )), and we will restrict 
ourselves to this case in this subsection. 

We first give a useful representation for E^'^* (a general representation for E'^'^* can also 
be obtained HI). From Eqs. ( glsD and (12^211) , we get 



27r2 



-P 



2iu 



bn,/ J {|^-p/2| |e + p/2|} 



l+2iu 



-P e 



l+2iv ■ 



K^' ^ {i^'i |e' + p|} 

Introducing the Feynman representation of the integrand 

1 _V{l + 2iv) f^lail-ay^'^-^da 





(2.24) 



r)l+2iu 

ot 



where 



L)„ = (1 - + + p)2 = + ctpp + (?a{\ - a) 

and setting R = B! + ap, Eq. (|2.24|) now reads 



(2.25) 



(2.26) 



+ p2a(l - a] 



T+2U7- (2-27) 



The integration with respect to the angle {q,R) leads to a Bessel function Jo^qR). The in- 
tegration with respect to = |^| can then be carried out, using the formula (6.565) of Ref. 
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J,{qx)x''+\ a^-^'q^'K.^^iaq) 
■dx — 



[x2 + a2]/^+i 



(2.28) 



This finally yields 



jq.p/2 (q 



2iv 



V2y r2(i + zz/)^o 



da[a{l — a)] 2e ^'^i'^K_2iu [qpy ot{l — a 



This last integration can be performed and yields (see appendix |A.6| ) 



(2.29) 



E'/ip) 



J — iiJ \ .6 1 J ~ iu \ .6 
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(2.30) 



where \l/ is the angle {q,p)- We will not need this explicit form in the following calculation. 



Let us now compute the contribution to the density of dipole ( [2.23| ) corresponding to 
n = 0. Since the longitudinal and the transverse degrees of freedom can be easily separated 
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(see formula (|2.23|) ), we will forget the exponential term depending on the rapidity Y, and will 
restore this dependence at the very end. We thus define n^n,u}{x,x',b) such that 



n{x,x!,Y,b)= J2 J Y-n{n,u}{x,x',b)exp(-^^x{n,i^)Y') 

In the asymptotic case we are interested in, 

, ~ , [ dv , f^aNc 
n{x, x,Y,b)r^ / — rijo.z.} {x, x , b) exp x(0, i^)Y 

J ATI \ TX 



(2.31) 



(2.32) 



We can now express the distribution of dipole, taking into account Eq. ( p.32|) . When y ^ oo, 
one can use a saddle approximation for the z/— integration. The function x(0, v) is maximum 
at = 0, which thus defines the saddle-point. We will thus expand ns^^ y^ around = 0. Using 
the expression (|2.29| ), ris^Q^y^ then reads 



n{o^y}{x,x' ,b) 



IGtt^ X 



1 



,o / X / N / . da C d(3\a{l-a)]~2\(3(l- (3)]-\ 

|6o,.Px'r2(i + zi/)r2(l-zz/) ^ (2vr)2 7o io ^ ^ ^'^ "^'^ 

X e -'2 " 2a)-iq.-{l - '^P) -iq.b^^^^^^^^^ _ a))K_^^,{qx' ^ ^{l - /?)). (2.33) 

Setting v_ = |(1 — 2a) — =-(1 — 2/3) + 6, the integration with respect to the angle {q,v) gives a 
Bessel function jQ{qv). The integration with respect to q can be carried out using the following 



formula 34 



/ x'^'^^ Kx{aix)Kx{a2x)JJa3x)dx 
Jo 



(^Y^ 'P77'{z){z' - l)-^^^r(/i + A + l)r(/. - A + 1), (2.34) 



2V 2 \a1a2) ^ 
X _j_ ^2 _j_ ^2 

where z = — and P is a Legendre function. Here = 0, A = 2iu and the 

2 aia2 
corresponding Legendre function is 



i 








.4] 




z + (z^ - 1)2 




z+(^z^ -1)2 





(2.35) 



In the case we are interested in, ai = xJa{l — a), a2 = x'Jp{l — /?) and 03 = b. In the domain 



6 3> X, x' (where the total cross-section gets its dominant contribution: see subsection 2^ 

52 1 



^-^[«(l-«Ml-/9)]"2 »1, 



(2.36) 



and thus 



J qK2iu {qx^a{l - a)) K^2iu [qx ^ (3{l - (3)^ Mqv) dq ~ ^^r(l + 2iu)T{l - 2iu) 







"62- 


62 ' 




xx' 



2iu 



[a(i-a)/5(l-/5)]— - 



xx' 



-2iu 



[a{l-a)(3{l-P)r 



(2.37) 



Using Ref. [|9 



7rV2 \h 



|2 



n,i/ 1 



z/2 + nV4 27r 

which gives for ^{o,;/} 



2 



f2.38) 



, , ^, 2zi^ X 1 r(l + 2zi/)r(l - 2ii/) ^ /-"^^^ 

(^— j [a(l - - /3)]-2- - j [«(! - - J . (2.39) 

The integration with respect to a and (3 can now be carried out. This yields 

Using the doubhng formula, one finally obtains 

This can now be inserted in Eq. ( |2.32| ) to calculate the large Y limit of the density, when 
}? jxx' ^ 1. We develop x(0, v) around i/ = 

X(0, 1/) ~ 2 ln2 - 7C(3)z/2 + (2.42) 
and take into account that x(0, —v) = x(0, u). One has then to compute 

?7,(x,x , r ,0) ~ -— / exp 21n2y 7n3)Y u + 2zz/ln . 

x' ¥ J -00 271 n \ n IT xx' ) 

(2.43) 

After making the replacement v' = v — % {Xtv^^^ j (^^^^7C(3)y'^ , the gaussian integration 
with respect to v' finally yields 

n(x, x', F, 6) ^ -7T3— ^ —hr^ exp In 2 F exp - ,Ln. LJ ■ (2-44) 



4^^-'(|«.iV.C(3)F)^/^ ^ - J I ^C(3)F 



This result is valid in the domain 



'''''''-7((3)Y«lu'-^«l. (2,45) 



TT XX' 



This expression corrects the expression (8) of Ref. |jT6|] when one considers the distribution 



in transverse space. Numerical simulations and approximate analytical calculations indeed 



confirm this factor 16 in the transverse distribution ||3^. Note that it is only in this asymptotic 
regime that n{x,x^,Y,b) has no angular dependence. 
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2.2 Exact result for the onium-onium cross-section at fixed impact 
parameter and equivalence with BFKL Pomeron 



Let us now compute the onium-onium scattering amphtude at fixed impact parameter according 
to the process displayed in figure |[ We denote x^i {xn) the transverse coordinate of the heavy 
quark (antiquark) making up the right moving onium and x^2 fhe coordinates of the 

corresponding quark (antiquark) making the left moving onium. These onia of transverse sizes 
2I1 = Xai ~ ^Ifei and X2 = Xa2 ~ illb2 scatter through the exchange of a pair of gluons between 
two elementary dipoles, respectively of transverse size x'^ and X2, located at b^i and 62 with 
respect to the reference point (which is arbitrary due to translation invariance). These two 
elementary dipoles are produced by the two heavy quark-antiquark pairs at a distance and 
62 from their center of mass. 




Figure 1: Onium-Onium scattering at leading order. 



In order to compute F^^^ as given by Eq. ( |2.7| ), we now use the result ( |2.15| , |2.16| ) for the 
dipole density and the following expansion of the dipole-dipole cross-section ctdd 



2al 



E 



n 



which is proved in appendix |A.2| . The full expression obtained is 

F'^'\x„X2,Y,b) = 



■IIM I 1,1 I 3i2 



X.2 



X / d w 



{2n) 



a2(16)2 g g g ^00 dUi dU2 



-«!=— 00 n2=— 00 n=— 00 



(2vr) 



(2vr) 



du 



(2.46) 
d'^Wi f d'^w_2 



1 2 



/ 2 

X2 



1 + (-1)" 



n-1 
2 



exp 



\ TT 



{x{ni,yi)Yi + x{n2,y2)y2, 
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CC-\ X 1 



^i.2^y--2,M2- Y--2y - ^y--2,-y--2 

X i?"'-* (fei + f - &i - f - I^) i?"-'^ (b,^t + &i + f - Ml, + 6i - f - I^Zi) (2.47) 

where we have used the fact that (Jj:,o{3li-,'^2ii^i ~ ^2) is translationally invariant and only 
depends on h^^^ = &2 ~ ^i- The quantum numbers rii, Vi and ?t,2, ^2 correspond respectively to 
the dipole distributions x'^, Yi, 61) and n(x2, a;2, F2, &2)- 

Note the existence of a term 1 + (—1)" which arises from the s ^ u symetry of the 
elementary dipole-dipole cross-section since the gluon is a massless vector boson (see appendix 
A.2 ). As a consequence only the even n will contribute to the scattering amplitude. 

The integration with respect to can be performed through the delta distribution. Then, 
the remaining expression can be strongly simplified when using the following orthonormalization 
condition for the functions E'^'^ (see Eq. (A. 16) of Ref. p9| ) 



/ 



p — E"'^(piO, P20) E"''^*{pio',p20') = anJn.mS{iy - fi) 5^(poO') 
IP12I 

+ (-l)"6„,,|poo'r'''*"(Poo'/Pm')X-m5(^^ + /^). (2.48) 



Note that this equation corrects Eq. (A. 16) of Ref. since the factor (—1)" in the second 
term of the right-hand side was missing. This term arises from the fact that 



-2^ - 1 \ ' , ... / 1 - -2^ ^ ' 



V.--i; <-^'"(— J . 

which has to be taken into account when performing the last transformation in Eq. (A. 19) of 
Ref. Applying the relation ( |2.48| ) for pi = fo^ + P2 ^ -i ~ ^' /^o ^ — 1 ^^'^ P-O' ^ — 

on one hand, for p'^ = &i + &j„t + ^) = ^1 + hint ~ Pq = 2^.2 + 1 and p^, = w on the 
other hand, and using the equality c?^&i rf^x'^ (i^Xg = dP p^d? p^cP p'-^d? p'^, the integration 
over these transverse variables gives 

17(1)^ n (16)2 +00+00+00 ^ ^ r,f,2 f ,2 



x/.v(..^4)( 



(27r)2 (27r)2 .f^^ ^-^ (2vr)3 J^^ (2^Y 

nl\ f . 1 + (-1)" 



' 4 i r ^ 4 i ^,2^f.-nn^„2^r..+i^2 



+ ( ^ I 1(1^" + 



X exp f ^^^^(x("'i, ^i)yi + X{n2, ^^2)5^2) 

\ TT 



/ ^ _ ^ \ ni 

ani,i/i5rM,n5(«^l " «^) ^^^(m " ^l) + ,1.1 " Ml T^"'^*''' — 7 _„5(z/i + u) (-1)"' 



— W 



1 / 

n 



an2,;.2^n2,n^('^2 - '^)5^(l£-m2 -^) +^n,i.|l£-I£2 -^1 ^ ^"^ \ 5 IT 1 "^"-"2 



X 5(i/2 + ly) (-1)" 



_ y;^, _^ _ ni,) E"^'-^ f ^ - w;2, - ^2') • (2.50) 
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The two terms involving the normahsation factors bni,ui and bn,i^ can be reexpressed using the 
fact that E'"''^ and E"'''^* are related by the following relation, which corrects Eq. (A. 12) of Ref. 
^ (see appendix |A.7|) 

h* / * \ " 

^'^'^*(Pio'P2o) = ^J ^Vo'i?"'^(p,o"P20')l/^oo'|-'^'^^ f (2.51) 



Note that this relation arises from the equivalence between the two corresponding representa- 
tions of SL{2,C) [0. The integration with respect to Wi and W2 can then be performed. It 
gives 



b* 



\ 2 


Ml, 


2 






K2 







)iM-ar«"'f4^)""(-i)- 

/ \w* — WiJ 

-w). (2.52) 



and 



-2 ^2 2 -V'- ^ 



i^-"'-^*(^y-2^ + &,-y-2^ + &j, (2.53) 

where we have used the fact that ii^~"'~^*(pio, P20) = -^"'''(PiO; P2o)- Since b_n,~u* = &n,y (see 
Eq. (p.22|) ) and a^n,-u = CLn,u (see Eq. ( p.38|) ), the contributions of the four terms obtained 
after expanding the brackets in Eq. ( 2.50|) are identical and one finally gets, using Y = Y1 + Y2, 



F«(a:i,X2,y,6) = -^(2vr)« j'^Ju j d?w[u' + 



+ (¥) -^ + (^) 



Z/2 



X exp ('^^x(^, i^)y] -w.-^- m] E'^'^ (^-w + b.-^-w + b). (2.54) 



TT 



Some comments are in order about this result, obtained without any approximation. First, 
it is clearly independent of the choice of the reference frame, since the result only depends on 
the total rapidity Y . In appendix [A.4| we explicitly show that this formula describes in an 
equivalent way, in the laboratory frame of the left-moving onium (which is defined to be the 
frame where this onium is moving relativistically, but has not enough rapidity to reveal its soft 
gluon contents), the scattering of this non evolved heavy quark-antiquark pair off one excited 
dipole at a distance b from the center of mass of the fast right-moving onium. 

Second, and more importantly, it explicitly proves the exact equivalence between the 
dipole and the BFKL approaches at leading order. Indeed, taking into account form factors 
when coupling the t-channel bound state of reggeized gluons to the external quark-antiquark 
pairs (see appendix |A.2| for details) and the difference of definition of amplitudes {Adipoie = 



^Abfkl), we should have the following relation between the dipole and the BFKL result 

F^^\xi,X2,Y,b) = J ^exp(a;f) [/^(x„i,Xbi,x„2,^fe2) + /^(^6i,^ai,^2,^b2)], (2.55) 
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where fujixai: 2lbiy iLa2j il±b2) IS defined by equation (26) of Ref. ||29| 

Juj[Xal, Xjji, Xa2y ^2) — / d SLo \ ^ + ~r ] ~/ ; . 2\ / 



X 



2asNc 



X{n, v\ 



(2.56) 



In this formula the integrand gets a factor (—1)" when permuting Xai ^ Xf,i. Performing the 
changes of variable w = Xq — -"i^-^i and (n, z/) (— n, — i/), one then recovers exactly the 
expansion ( [^.54D with the factor 1 + (—1)", which proves the result. Note that the equivalence 
between the BFKL and dipole kernel can be also proven by comparing the real and virtual 
graphs in covariant and light-cone quantization. The result is that the sum of real and virtual 
contributions is identical in both case, although each of these terms differs. Thus, this result is 
true only for inclusive quantities |1T7 |. 



Defining F^l^^^^ as 



F«(xi,X2,y,6)= E /_7^^{nUfe,2l2,&)exp(^x(^,z.)f 



/2a, AT, 



(2.57) 



and using Eqs. (|2.19|) and ( p.31| ), equation ( |2.54| ) can be rewritten as 

.(1) . ,^ T^aW?. i + (-r" 



n+l 
2 



(2.58) 



In the asymptotic regime where one can keep only the term corresponding to n = 0, this 
relation simplifies to 



^{0,1} (^1,^2,^) 



2^2 



s-^2 



n 



{0,u}{Xi,X2,b). 



(2.59) 



Computing F^^\xi,X2,Y,b) by a saddle point method as we did for n{xi, X2,Y ,b) (see Eq. 
( |2.43| )), one has to expand the pref actor around = 0, which yields 



F^^\xi,X2,Y,b) ^ -4:iTa1xln{xi,X2,Y,b) 



,XiX2 ln(16 6^/xiX2) 

~ \ 3/2 



62 



(|a,iV,C(3)f 



exp 



IT 



\n2Y\ exp 



in the domain 



^^7C(3)F < In < 1. 

7T X1X2 



\n^{lQb'^/xiX2)\ 
^C(3)r j 

(2.60) 
(2.61) 



This result, which differs from Eq. (10) of Ref. ||T6[ by a factor 16, is in aggreement with 
numerical simulations ||35[| . 
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2.3 Calculation of the onium-onium total cross-section 

In this subsection we compute the onium-onium total cross-section. It is related to the onium- 
onium cross-section at fixed impact parameter by 



(2.62) 



We show that, provided the elementary dipole-dipole cross section is integrated over distances 
but not averaged over angles, one can get the onium-onium total cross-section by a simple 
direct calculation. 

Combining formulae (|2.7| ) and ( p. 62 ), one obtains 



Ftot {^Lii !!!L2i ^) 



2'Kx'i 2^X2 

^^72 o ~% n{xi,^^, Yi) n(x2, X2, Y2) I d^{U2 - Ui) ^2' ^'1 - ^2)' 

(2.63) 



where n(xj,x-,Fi) is the integrated density of dipoles 

n(x, x', y) = j d"^!} n(x, x/ ,Y ,b). 
Thus, since b and q are Fourier conjugated. 



(2.64) 



+00 



n=— 00 

+00 



du 



\x\ 
\x' 



X(n, u)Y 



dv 


\x\ 


/ x*x'^ 


n/2 

) 


x' 


-2iu 

exp ( 


2^ 


\x'\ 


\xx'* J 




X 



TT 



-x{n,u)Y 



(2.65) 



where we have used Eq. (|2.23| ) and the expansion ( |A.32| ) of E'^'" for g 0. The integration 
of (Jdd{2±i,2±2^1^i "^2) with respect to the distance 62 — 6']^ is performed in appendix |A.2| and is 
given by formula ( A.33|) , which still depends on the orientation of the elementary dipoles. One 
now gets for Ff 



(1) 

tot 



Ftot {SLi^ 2l2i^) — \ 51 



X 



Xi I X\Xi 



Xi \XiXi 



, +2? r+^dui +2? r+^du2 f d^x[ d^x'2 

^ y 00 27 ^ / ~ 

ni=_oo"'-°° n2=-oo 
2iu2 



00 2t[ J 2nx'i 2ttx2 



l2 



^n.i/2 




— 2wi 


X2X2 


n2/2 


X2 




Xi 


X2 


\ 3^2 3^2 / 




X2 



exp 



2aNr 



X a 



2 ■^1-^2 



' I / * ^ 
^1^2 



TT 

-7i/2 / 2iv 



+ (r^X] ^2 



{x{ni, Vi)Yi + x(n2, 2/2)^2) 

(2.66) 



Integrating with respect to x'^ and x'2 leads to a lot of 5 distributions since 



]_ r dx[dx[ !n_ii_i(^,_^) _ 

27ry 2x[x[* ' -^dn„nd[u, u) 



(2.67) 
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and a similar integral over x^- Defining F^^^^^^^-^[x_x,x_2) and X2) via 



and 



vr 



?^{n,i.}(^i,X2) = 2^ / — n{„,i.}fe,X2) exp ( — —x{n,u)Y 



this finally yields 

-^tot{n,i/}fe'^2) 



TT a t 



i + i-r 



TT 



(2.68) 



(2.69) 



4 ,^2 



2 ^2 



s -^2 



l + (-l)" 



T Fll F2I 

XIX2 / 



(2.70) 



Comparison with Eq. ( p.58|) provides a check of this calculation when integrating both sides 
with respect to b. In the asymptotic regime, corresponding to n = 0, 



2^2 



natx 



n{o^y^{xi,X2), 



(2.71) 



which could be obtained from Eq. ( |2.59| ). Using Eq. ( p.65| ), the corresponding asymptotic 
integrated dipole distribution reads 



ri(xi,X2, r) = / — exp x[n,y)Y 

2nX2 \X2J \ TT 



(2.72) 



which gives, in the saddle point approximation at large Y , 



1 xi 1 



In 2 y 

23^2 y'la,iV,C(3)F 
Thus, one gets for the total cross-section 

F^^\xi^X2^Y) ~ —ATia1x\n{xi^X2,Y) 



exp 



ln^(a:i/a:2) \ 
^C(3)rJ 



(2.73) 



exp 



—2tx a Xi X2- 



ln2r 



, j ln^(xi/x2) I 

|«.iV,C(3)F """^l i^C(3)rJ 



(2.74) 



in agreement with formula (26) of Ref. [T^ . 

Let us integrate the scattering amplitude F^^^ (6) with respect to the impact parameter b 
in the domain (|2.61|) where formula ( p.60|) is valid (neglecting the fact that the upper bound is 
not infinite). The corresponding integration then gives 



F^^\x^,X2,Y) ~ -27T 



exp 



«, Xi X2- 



\n2Y 



IT 



asN,C{3)Y 



\ ln^(xi/x2) 1 



(2.75) 
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which is identical to Eq. ( p.74|) . From the gaussian distribution obtained in Eq. (|2.75|) , it 
is clear, comparing with Eq. ( p.60| ), that the total cross section at BFKL order is dominated 
by impact parameter configuration much larger than the transverse sizes of the two scattering 
onia, corresponding to 

■ 1662 \ lUa^N, 



In 



X1X2 



-cm. 



(2.76) 



Note that this dominant contribution is inside the domain (|2.61 
tions are much more central than what was claimed in Ref. 



These dominant configura- 
It confirms previous numerical 
simulations [^. Thus, the calculation, based on perturbative QCD, is expected to remain valid 



m 



for high values of Y. 



3 Electron-Onium Deep Inelastic Scattering 

In this section, we perform an analysis of — onium deep inelastic scattering at low Xbj, based 
on fcy-factorization and dipole color model, as illustrated in figure 0. Our aim is to compute 




Figure 2: /cT^-factorization and dipole model applied to e — onium deep inelastic scattering. 

various structure functions in the small Xf,j regime, where the onium wave function is dominated 
by a color dipole cascade. 



In the Regge limit, one can apply the fc-r-factorization tool |2y, ^ ^ in order to extract a 
photon of virtuality off an onium. It involves the elementary Born cross-section a^g/Q'^ of 
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the process 7 g{k) qq. Here the gluon is off-shell, quasi transverse, with a virtuality k"^ ~ fc^. 
One also has to introduce the unintegrated gluon distribution density at a factorization scale 
(5oj which is related to the usual gluon distribution by 

G{x,,,Q\Ql)= d^kJ^ix,„k,Ql). (3.77) 



We first deal with a dipole of transverse size xqi, which can be either part of a heavy onium (i.e. 
a heavy qq pair) or extracted from a proton as will be emphazied later. The /cT-factorization 
implies, for the total 7* — dipole cross- section cr^* , 

Q^a^,ix,,,Q^;xl,) = J d^k £ ^ a^g{x,,/z,kyQ^) J^{z,k;xl,), (3.78) 

J-'{z,k;xQi) being the Fourier transform of J-'{z,k;QQ) in transverse space. We next evaluate 
by coupling the gluon to the softest dipole which arises in the cascade. This is achieved by 
using a second fcr-factorization. It involves the elementary Born cross-section Cgd/k"^ of the 
process d{x) g{k) — ^ d{x) for a dipole of transverse size x and a soft gluon of virtuality fc^. 
This /c^-factorization can be expressed by 

^J^{z,k:,xl^) = I ^ 2 / —r^^^oi^^^'^-l)^9d{z/z\£^)5{z/z' -I), (3.79) 

J yZTl J X J z Z Z 

where the distribution density n{xQi,x,z) was defined in section As previously zip+ is the 



light-cone momentum of the quark part of the dipole (xqi). is computed in appendix [A. 3 
using eikonal techniques. Defining Y = \nzi/z, one gets for the dipole-photon cross-section 



Q'^a^,{xbj,Q^;xl^) = J d^kj^ ^igi^bj/ z,^) J ^n{xoi,x,Y) 

Let us now compute the convolutions in longitudinal and transverse spaces. As in section |], 
we introduce a double Mellin-transform in these both variables, namely 

nixQ^,x,Y) = J ^e'^^n^{xQ^,x) (3.81) 



and 



{xo„x) = (3.82) 



Here the Mellin variable in transverse space is 7 = | -|- zz/, which is introduced here rather than 
v since it plays the role of an anomalous dimension. We consider only the dominant Regge 
trajectory, that is n = 0. Formula (|2.17|) then reads 



where 



X(7) = X (0, ^ + ^y) = ^(1) - i^(7) - 1^(1 - 7) (3.84) 
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The quantity has been calculated for different polarizations of the incoming photon in . 
We introduce the corresponding double Mellin transform of this cross-section, namely 



or equivalently 



4vr a.../..(7) = 1^ I ) I ) (3.85) 



= 4vr'«e.^ f^(^] ' -^-(7), (3.86) 



with 



The expression for cr^jj, now reads, after performing the convolution in longitudinal space, 

a^Nr, A ,0. f d'-y' f d'-y f du f , z-i \ / XqiX"^^ 



J J ztn J 2i7T J 2m y x^j J 



(4 o.^.[x,,M ,a;oij-47r «e.™^y -^-y 2z7r 7 2^7 2^7r^'^"\" a;, J V x ) 



C^-^X(7) Y i U^j ^ ^ 

The integration with respect to the polar angle of x leads to a Bessel function. One has then 
to integrate over x, namely 

/ ^47r(l - Mkx)) — = Anikxoif'' , ^ = 47r(fcxoi)%(7). (3.89) 

J X \Xoi/ 7 1(1 + 7) 

The integration over k gives 7 = 7'. Since ^^(7) (formula ( ^.83[ )) exhibits a pole at Up = 
^^^xil)^ the a; integral finally yields 

J^<i.M) ^ ^/|^/.,W^(4.Q=)^exp (H^XW..^). (3.90, 

In the regime we are interested in, cjp <^ 7 and the dependence of h^^ (7) on Up can be neglected, 
replacing Up by |20|. We then get for the total cross-section Y ~dipole{xQi) and for the related 
structure function (we assume R = Fi/ Ft to be small in the relation between the total cross- 
section and the structure function) 



(3.91) 

The dipole (xqi) being part of a bound state, for example extracted from an onium, one has 
now to average with respect to the coresponding wave function, $*^°)(xoi, zi). The initial dipole 
state is supposed to be well localized in transverse space. Namelly, its transverse size, given by 
the scale M ^ which is defined via 

{M'r^ = J d'x,,{xl,ydz^<!>^'\x,„z,), (3.92) 
is assumed to be perturbative. One then obtains for the onium structure function 

(.,.Q-M=) = 2^ f^l^xp f?^,(7)l.^) . (3.93) 
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It is then possible to apply formula (|3.93|) to peculiar structure functions, namely 





(2a,N, 1 
exp xiV In — 

V Xbj 




(3.94) 



where -Ft(l) is the structure function corresponding to transverse (longitudinal) photons and 
Fg the gluon structure function. The coefficient functions 




as (r(i-7)r(i + 7))^ 

37r7 r(2 - 27)r(2 + 27) 1 



7 



;i+7)(i-i) 

7(1-7) 



(3.95) 



were computed in ref |20]. When is small, the 7— integration can be performed by the 



2' 



steepest-descent method. The corresponding asymptotic saddle point is located at 7 
which defines the BFKL anomalous dimension. Expanding the x function around |, we obtain 
a saddle point at 



7. 




where a 



TT 



(3.96) 



The approximation of expanding xil) around ^ is valid when 



a In 



Q_ 

M 



In^/ln— «1, 

M' Xbi 



(3.97) 



that is the small x^j, moderate Q/M kinematical domain. This yields 



Ft + Fi 



M 



(«P-l)ln--fln^^ 



(3.98) 



where ap is defined by Eq. (p,.2D. Thus, F2 depends only on 3 parameters, C, M and ap. 
Suppose we can fit F2 with this form. Then, get a prediction for Fq and R = F^/ Ft without 
any free parameter. Namely, 



Fg 
F2 



flT + hL 



37r7s 



1 



Ms 



r(2-27,)r(2 + 27,) 



and 



7=7s 



R 



h, 



a. 



i + |7.-|7Mr(i-7s)r(i + 7.))= 



7.(1-7.) 



hT^^'-'^ (1+^^)(1_^)^ 



(3.99) 



(3.100) 



where 7^ is given by the expression (|3.96| ) . Note that the overall non-perturbative normalization 
C does not enter R and Fg/F2. It is possible to apply this analysis to the proton This 
requires some additional asumptions when considering the coupling of the dipole cascade to the 
proton. It leads to a successful description of the HERA data It also provides a prediction 
for the gluon density based on the BFKL dynamics, and a prediction for the ratio R, using 
formulas (|3.99|) and (|3.100|) . All the previous annalysis was done by performing an expansion 
around the small Xbj behaviour of the BFKL Pomeron, that is by considering In 1/ Xbj as a big 
parameter, which leads to a saddle point close to 7 = 1 (see Eq. ( p.96|) ). A different analysis 
of the modification, due to the BFKL dynamics, of the anomalous dimension of the double 
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logarithmic approximation common to DGLAP and BFKL can be done by considering now 
InQ^/M^ as a big parameter. The corresponding saddle point of Eq. ( |3.94| ) is now around 
7 = 0. In Mellin space, the double logarithmic expression of the anomalous dimension is given 
by 7aj(as) = 3as/ (tto;), and expansion of the x function ( p.84| ) around 7 = leads to corrections 
given by powers of 3ac(/(7ru;). This method provides an extension of the domain of applicability 

As it has been seen in this section, the dipole model can be safely applied when the 
two scales of the process are both perturbative, as it is the case for — onium scattering. 
The application to — p scattering requires some assumptions for the coupling to the proton. 
Because of the well-known diffusion in transverse momentum space, such an application of the 
dipole model, although successful |2^, |2^, cannot be considered as a clean test of high- 
energy perturbative regge dynamics. A possible test of such dynamics could be based on single 
jets events in DIS or double jets events in hadron-hadron collision [Q. Another interesting 
test of BFKL dynamics would be the 7* — 7* events in — e~ colliders at high energy in 
the center of mass of the virtual photon pair and with high (perturbative) photons virtualities. 
This has been already proposed in the framework of the original BFKL equation ||3^. Such a 
process can equivalently be described in the dipole picture of BFKL dynamics. This will be 



developped elsewhere [40 



4 Conclusion 

In this article we have shown the exact equivalence between BFKL and dipole approaches 
for the onium-onium cross-section at fixed impact parameter. This proof relies on conformal 
properties of the dipole cascade and of the elementary dipole-dipole cross section. We have also 
obtained asymptotic expressions for the dipole distribution inside an onium and for the onium- 
onium cross-section at fixed impact parameter. These results agrees with previous numerical 
simulations. We also apply the dipole model to onium-e^ deep inelastic scattering, using the kx- 
factorization, and obtain predictions for various structure functions in the BFKL dynamics. The 
different elementary cross-sections used in this paper are computed using eikonal techniques. 

Relying on the same conformal properties, it should be possible to get analytical ex- 
pressions for the multipomeronic contributions to the onium-onium cross-section, which are 
expected to be important for large rapidities. 

From a phenomenological point of vue, the dipole framework could be applied to other 
inclusive processes. The application of this technique for exclusive quantities remains however 
an open question, due to the use of light-cone quantization, in which the intermediate states 
are unphysical. 
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A Appendices 



A.l Calculation of aDD using eikonal methods 

In this appendix we compute the dipole-dipole cross-section using eikonal techniques. In QCD, 
the eikonal current due to a fast quark of momentum p, responsible for the emission of a soft 
gluon of momentum k {k <^ p) and color a reads 

f{k) = -rgT^^--. (A.l) 
p.k + le 

Let us consider the scattering of two dipoles. Let Xq (x^) be the transverse coordinate of 
the free quark (antiquark) making the right moving dipole and (x'^) the coordinates of the 
corresponding quark (antiquark) making the left moving dipole. In this mixed representation 
where the fast radiating particule is represented in impact parameter space, and the radiated 
gluon is represented in momentum space, the eikonal current corresponding to a quark of 
transverse coordinate x then reads 

r{k) = -rgT^-^e-'3i-k (A.2) 
rj.k + le 

where t) = (l,0,0). This current is responsible for a term —ii°'"{k)ey when computing the 
amplitude of emission of a gluon by the quark. This terms arises when evaluating the evolution 
operator T exp —i J d'^xj°''^{x)Aau{x) due to the usual hamiltonian / (fix j"'''{x)Aau{x). Note that 
this current is completely described by the transverse size of the dipole and does not require 
any additional information about its internal structure. 
The eikonal current corresponding to the right moving dipole then reads 

j-(A;) = _z^T«_Jl_(e~^2Io-i _ e-'^^-^), (A.3) 

and a corresponding formula for the left moving dipole, replacing f/by?7 = (0,l,0). 

Let us now compute the graph Ai represented in figure ^. It reads, in Feynman covariant gauge. 




1 



P 



kl Ik 



.7 



1' 

0' 



p 



Figure 3: Contribution to the dipole-dipole scattering. 
^1 1 r d'^k 



X (-,/+(A:))(-./-(-A:))y£. (A.4) 
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The terms 2^2^^ related to the normahsation of the initial dipole states and ^'''^ and 

^^^^ are due to the fact that the considered amplitude is computed for amputed propagators. 
The color factor reads 

1 AT^ — 1 

Y.TrT'^T'TrT'^T' = -f-—, (A.5) 

c ab c 

which equals | in the large Nc limit. 
Using the expression ( [A.3| ), this yields 



f dH 1 1 



X e 



(27r)4 (A;2)2 k- + ie k+ - ie 
ixQ.k _ ^ixi-k^ ^Q-^ixo-k. _ ^-ixi-k--^ ^^-ixQ-h _ ^-ixi-h^ ^^ixQ-k _ ^ix^.k--^ g-^ 



One has also to consider the crossed diagram A2. In order to integrate with respect to k~ and 
k~^ we consider two equivalent representations of Ai and A2, obtained by changing the sign of 
k (see figure |[ Forgetting the phase factor for a moment, the forward amplitude A then reads 



»k kl 





A2 A2 
Figure 4: Contributions to the elementary dipole-dipole scattering. 



A = Ai + A2 



9^ 



(Pk dk- dk+ 1 1 / 1 



1 



4 7 (27r)2 27r 27r (F)2 2VA;- + ie -A;-+ie/ \k+ + ie -k 



1 



2e 



(A.7) 

The two last terms reads {—2TTi6{k )){—2TTi6{k^)), which finally yield, performing the integra- 
tion with respect to k^ and k^ 

■•2 d^k 



2 J (^2^2 



2-e 



2 g'fk-ikoi g~'£L-^l 



which was obtained in Ref. ||T5| by computing elementary Feynman diagrams. If we neglect the 
dependence of $(0:01,2;) and n{xoi,x,Y) with respect to the dipoles orientation (see appendix 
[A.5[ ), it is possible to average with respect to the angle of these dipoles. A then reads 

POO rjh 

A = -A7ialj^ j,(fca;oi)][l- Jo(Ki)] (A.9) 
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Using the identity 



1 + ln^ 



(A.IO) 



where a;< = Min{xoi,XQi) and x> = Max{xoi,XQi), this finally yields for the corresponding 
cross-section 



0-DD{xoi,x'(i 



01 J 



-2A = 2Tialx\ 



X- 



(A.ll) 



This average result is however sufficient when computing the total cross-section (see appendix 
Another representation of this elementary cross-section is very useful. Consider 

2 



2 J-oo [v 



2 , U2-^l -^2 



(A.12) 



When xi > X2 {xi < X2) one can close the integration contour around +ioo (— ioo) , so that 
one pick up the (simple) pole at u = i/2 {y = —ij^) . Thus, 



2 du 



a e ^< 



(A.13) 



where x< = Min{xi,X2) et x> = Max{xi,X2)- This finally yields the expected result: 



aDD{x[,X2) 



a: 



2 (z/^ + 



^^,{x[r'-{x'2y-'-. 



Note that this can be written equivalently as (see the following section) 



A. 2 Calculation of aDD{x, sf, b — 6') 



(A.14) 



(A.15) 



In this appendix we compute the elementary dipole-dipole cross-section at fixed impact para- 
meter. As in the appendix A.l , we consider two dipoles of transverse sizes x = Xj^ — x^ and 



—0 



,1 xU+x- 



x! = x!i — Xq, whose centers are situated at 6 = i^lli and U — ^ . 

Let us compute the non-forward scattering amplitude of these two dipoles. In the high- 
energy limit —t<<s = 2p.p', the exchanged momentum is quasi-transverse 



g- (0,0,g) 



(A.16) 



The graphs to be computed are displayed in figure |^. These graphs are computed using the 
rules of appendix |A.1| . Since q~ = = 0, the expression for the corresponding amplitude is 
similar to the one obtained in equation (|A.8| ), except for phase factors. It reads 



A = 

2 



r d'^q 


f d'^k 1 


1 {2nyJ 





^ixQ.k _ ^ix^.k\ L~ixQ.{k + q) _ ^-ix^-ik + q) 



-IX r\ . k 



(A. 17) 
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kl fk+q k+qf ik 



Ax 





A2 A2 

Figure 5: Contributions to the elementary dipole-dipole scattering with non zero exchanged 
momentum. 

Note that this expression is overall translationally invariant. However, it will depend on the rel- 
ative distance between the two dipoles and on their angles. Expanding the previous expression 
and integrating with respect to the angles of k and k + q, one obtains several Bessel functions 



A 



where 



•y {J(|Xo-Xo/|) -/(|xo -Xo/|) - Xr\}^ , 



r+oo 

H^) = J Jo{kx). 



(A.18) 



(A.19) 



computed (see Ref. fl^ ) 



p is an infrared cut-off which regularize the divergency at A; = 0. In the limit p — > 0, / can be 

ipil) + In 2 — In a; — In p. 



/ = lim 



k 



— Jo{kx) 



k 



(A.20) 



When evaluating A, the constant ipil) + ln2 and the infrared divergent term Inp cancels. A 
then reads 

' 2 



A 



In 



Defining 

^ _ Xq + Xi 

2 

the non-forward cross-section finally reads 



^0 


2i0' 




1 ^ 2 
^I'l 1 




—1' 1 


1—1 


3iO' J 


and 




_ 2i0 


+ ^1' 






2 



(^DD{x,x^,b — U) = — = 77 



(A.21) 
(A.22) 

(A.23) 



as quoted in Ref. Note that it only depends on the relative distance h — U because 

of the overall translational invariance. This expression can also be obtained by evaluating 
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this scattering amplitude in the laboratory frame of one of the two onia. In this frame, this 
amplitude is expressed as an eikonal phase (computed in terms of a Wilson- loop), due to the 
change of the wave function of the slow moving onium in the color field of the fast moving onium 



26|. Let us now show that this cross-section can be expanded on the basis of the functions 



E""'^. In the standard Regge calculation, instead of considering the scattering of two dipoles, 
one considers the scattering of two gluons of momenta —k and k + q (see figure |A.24| ). In the 




Figure 6: Gluon-Gluon diffusion in the Regge framework, 
lowest order approximation, the contribution of this graph is 

s\k-k!) 



A 
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k\k + qY' 



which gives, after Fourier transform in impact parameter space |¥T 

a(^o,^i,^o',^i') = (27r)^ln(|xoo'|A) ln(|xn'|A), 



(A.24) 



(A.25) 



where A is the mass of the gluon introduced in order to remove the infrared divergency. Since 
one is interested in the coupling with color neutral states, this expression can equivalently be 
replaced by 



3Ioo'2ill' 



3Ioi'2iO'l 



In 



3Ioo'2ill' 



3Ioi2iO'l' 



(A.26) 



because of the conservation of color current [^. In the case of the dipole model, instead of 
labehng the gluons, one labels the dipoles. Thus, the contributions of the figure ^ should be 
equal to the contribution of figure |^ plus permutations (xg ^ ^i) and (xq' ^ Xr) (except for 
normalisation factors). Indeed, one can check that 



crDD{x,x!,b- ^) 



(27r)2 

2aJ 

'^(^O' — 15 — 0'' — 1') ~^ (—0 ^l)]' 



(A.27) 
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Using the expression (26) of Ref. in the Born approximation, that is making g = 0, one 
finally obtains 

i2\ 



2a^ +°o 
(Tdd{x, x, h-h!) = — 



T 



1 + -If 



,2 I / n+1 



(A.28) 



smce 



E'^^'ixho, Xao) = i-irE^'''{xao, x,o). (A.29) 

One can check on this expression that when integrating over U — h and averaging over 
angles, one recovers the total cross-section ( |A.14| ). Indeed, using the mixed representation 
(|2.21| ), one obtains 



1 + (-1)" 



(2vr)2 

X e^'''*{x')e;;'''{x). 



+ (^) 

(A.30) 



Thus, 



anoix^x') = J d^Ql — h) crDD(x, x ,b — l/) 



I +00 

' 16 ^f-L 



dv 



^ lim 



(A.3i; 



Using the expansion (A.l) of Ref. [^, one obtains 
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"/2 ^ -2i;. 
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1 1 — 4ii^g— i(5(rt,i') 
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1 1 —iiu^—i5{n,v) 



(A.32) 

where e*''^"''^^ is a phase given by Eq. (A. 2) of Ref. [^. When g ^ 0, only the two first 
terms remain. Due to the symetry of — ^Xi')^ — : nrv when v —v or n ^ — n, the 



contribution of these two terms is the same when summing over n and integrating over v. Thus 



aDD{x,x) = WJ-^-^ 



16 



dv 



p2 (n+lY\ yX'X 





n/2 
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/ xx* ^ 
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\ x*x' , 
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-2iu 



Averaging over the angles, the only remaining term is n = 0. Thus, 

n/2 r+co 1 



aDD[X,X 



l+2iu \ I \l-2iu 



c/z/|x|™|x' 



2 ' 



(A.33) 
(A.34) 



which is identical to Eq. ([A. 131). 
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A. 3 Calculation of cFgd using eikonal methods 

In this section we compute the elementary Born cross-section a^d/k"^ of the process 

d{x) g{k) d{x) (A.35) 

for a dipole of transverse size x and a soft gluon of virtuahty h^, in hght-cone gauge. This process 
is illustrated in figure Consider a dipole of transverse The corresponding 



k, A, a 



k, A, a 



+ 



Figure 7: Amplitudes contributing to the elementary dipole-gluon cross-section. 



eikonal current has the same expression than in covariant gauge ( |A.3| ). This would not be the 
case for a more complicate system for which the current would be rotated in color space (see for 
example Ref.|^). In light-cone perturbation theory, fc^ = 2k~^k~ — k^ = 0. The corresponding 
current then reads 



(A.36) 



The Born cross-section then reads, summing over color and polarization of the emitted gluon 
and averaging over the color of the dipole. 



e 2(27r)2Ar, 



Tr\j^t 



-\\2 



^ Y,g^Tr{T^T^) [e'^^-^ - e''^'^^ (^(e^^-^a _ e^^-^f) 4 (^) (A.37) 



c A=l,2 a 



Here one takes into account the fact that the dominant polarization of the soft emitted gluon 
in light-cone gauge corresponds to 



(5, 



k. 



(0,^,0) 



in the soft limit. This finally yields 



k' 



2(27r)^ 



Jk.x 



-ik.x 



asNc 
7r 



Jk.x 



-ik.x 



J2- 



(A.38) 



(A.39) 
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Let us make two comments. First, the dominant polarization of the exchanged gluon is lon- 
gitudinal, due to the expression ( |A.38| ). However, the polarization tensor in light-cone gauge 



reads 

d^^ = g''" - ^ V • (A.40) 
T].k 

In the Regge limit, since ji, is proportionnal to fj, the relevant component of this tensor is d'^~ , 
and since e"*" = in light cone gauge, it turns out that e^d^' = e*(i*~, and thus the exchanged 
gluon has a physical polarization, as it is expected in high-energy factorization. 
Secondly, when integrating this cross-section over k, one can recover the dipole emission kernel 
which was constructed in Ref. 



A. 4 Calculation of the onium-onium cross section in the laboratory 
frame of one onium 

In this appendix we illustrate the frame invariance of the onium-onium cross section by an 
explicit calculation in the laboratory frame where the left-moving heavy quark- ant iquark pair 



is close to be at rest. Let us first consider as in section |2.2| the scattering of two heavy quark- 
antiquark pairs of transverse sizes Xi and X2 which scatter through the exchange of a pair of 
gluons between two elementary excited dipoles of transverse sizes x'^ and x'^. The corresponding 
expression for the scattering amplitude F^^^ is displayed in formula ( ^.47] ). Using the same trick 
which led us to formula (|2.50|) , we integrate over d'^b^.^f-d'^af^ and d'^W2, which yields 

F^'\x„X2,Y,b) = -2^^- EE/ 77A-3 du I 



27TX[ 



l2 




{2n) 



00 J Xi 



' 2^n^\ l + (-ir f 2asN, ^ , ~ ~, 

^ + -r] 7 : — r^w : — ^xp [x[ni, yi)Yi + '^)>2. 



^1 



TT 



2 2 V2 2 



Ml 



X E-''^* l^b^ + ^-nj,b,-f~wj E-''^ ^^-nj + b,^-w + bj (A.41) 

Integrating over d'^gcf^d'^bi and d'^w would give the result ( |2.54| ). It would show in particular 
that in the previous formula rii can be replaced by n and z/i by u. Since Y = 1^1+1^2, the 
previous formula can be written as 

F^'\x„X2,Y,b) = -^l d%l £|^n(xi,x'i,l^,&i)^DD(2l'i,ai2, (A.42) 

which is exactly what one would obtain when computing this process in the boosted frame where 
the left-moving onium is close to be at rest. Indeed, let us consider the dipole distribution inside 
the left-moving onium in its laboratory frame, namely n{x2,x!2,Y = 0,62)- Using expression 
(l2l9| ), it reads, using the completeness condition for the functions E'^''^ (see equation (25) of 
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Ref. dH)) 
n{x2,x!2,Y = 0,62 




(27r)3-^^ "V 2 ; " ^ ' 2 

= 27^x1 6\b2) 6\x:2 - X2). (A.43) 

Calculating the scattering amplitude F^^^ with formula ( |2.7| ) in this laboratory frame, one gets 



X2, 6, y) = -I / c;^&2 d\y2 - &;) - 62 - U, + &'2 - &) 



X n 



fe, ^i, F, 61) 27ra;2 ^^fe) ^^(^2 - ^2) crDD(^i, ^2, - U2] 



= -\j d\ J ^^n(xi,x;,y,6i)az)Z5(x;,X2,6i -6), (A.44) 

where we have performed in the last step the integration with respect to (i^X2 rf^&2 '^^(^2 ~ ^'i)- 
This result is identical to Eq. ( |A.42| ). 

A. 5 Approximate calculation of F^^^ in the asymptotic regime 

In this appendix we calculate the onium-onium scattering amplitude at fixed impact parame- 
ter in the asymptotic regime. In that case the distributions of dipoles n{xi,x[^,Yi,hj) can be 
approximated by keeping only the contribution = to formula ( |2.23| ). We show that even 
in this limit it is important to consider the exact expression ( |A.23| ) of the elementary dipole- 
dipole cross-section. In Ref. []T3[, the onium-onium cross-section at fixed impact parameter 



was computed considering the interaction of two elementary dipoles of transverse size x'^ and 
x'l, situated at the same point 6'^ = 62; interacting through the elementary total cross-section 
udd^x'^^x^), and the dipole distributions were averaged over angles. Indeed, the elementary 
cross-section at fixed impact parameter cr£)£)(a;'i, Xg, fe'^ — U2) is dominated by configurations 
Ui ~ ^25 since c'"dd(^15 aig, 6'^ —Ih) ~ \y \i when —^2! is large (this can be seen by expand- 
ing formula ( [A.23| )). Thus, it is possible in first approximation to make 6'^ = 62 and to replace 
Ui ~ U2) by its integral over 6'^ — 62? which is equal to ct£)£)(x']^, Xj) when averaging 
over angles (see appendix |A.2|) . Note that in the asymptotic regime where it is possible to keep 
only the contributions ni = ?t,2 = in the dipole distributions, the elementary cross-section will 
only contribute through the quantum number n = 0, since this quantum number is conserved 
(see Eqs. ( p.47[ - pl54|) for example), and thus the averaging over angles in the elementary cross- 
section is automatically implied (see Eqs. ( A.33 - AT51 )). However, we will see below that the 



approximation of integrating the elementary dipole-dipole cross section over impact parameter 
is exact only when calculating the total cross-section. 

We thus calculate first the cross-section at fixed impact parameter in these approximations 
and compare the result to the exact formula ( |2.59| ). Let us consider two heavy quark- ant iquark 
pairs of transverse size and X2, which scatter through the exchange of a pair of gluons 
between two elementary dipoles, respectively of transverse size x'^ and X2, and at a distance 
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from the center of the quark- ant iquark pair bi and 62- In this approximation, the corresponding 
scattering amphtude F^^^ given in Eq. ( p.7| ) simphfies in 

~ yoo (^j^' c/x' ~ 

F^^\xi,X2,Y,b) = I — r^— ^rf^&iri(xi,x'i,6i,Yi)n(x2,X2,6-6i,F2)o^z)Z)(a;'i,X2). (A.45) 



The rapidity and Y2 are such that Y = Yi + Y2. F^^^ involves the elementary dipole-dipole 
total cross-section, which has been evaluated in |]T3[, and which is calculated in appendix [A.l| 
using eikonal techniques. For two dipoles of transverse size x'^ and Xg, the elementary forward 
dipole-dipole cross-section reads (see Eq. (|A.11|)) 



2 ™2 
< 



X2) = 2na1x 

where x< = Min{x'i,X2) et x> = Max{x[,X2). 
We write F^^'' given by Eq. (|A.45|) as 



1 + ln 



(A.46) 



/CXD 7 / 

— F«(xi,X2,6). (A.47) 
-00 ZTT 

Using the expression (|A.14|) of the elementary dipole-dipole cross section ann, Fj)^^ then reads 



2 /^2+l 



21 — / — exp ( ^— (x(0, Ui)Yi + x(0, z/2)>^) 



—00 ' —00 Zjc 



TT 



27rxf J 27rx2^ 



d\i / (i^62^^(^-k-^)'^i'i(2;i,2;i,6i)n^2(x2,x2,62) 



Using the representation 



— exp(-ig.(6 - 61 - 62), 
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(A.48) 



(A.49) 



and the expression ( p.23|) of n(x, x', one gets 
^(1) _ Tra 



^ 1 f (Pq dvi dv2 f2af,N,, 
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27rxx^ J 211X2 
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(i x'^ f d X2 ^—iq.h ^ ^0,1^1 

Xi 



>l+2w / / \ l-2ii^ 



(A.50) 



Let us compute the integrations with respect to x^. Defining 

where ei = 1 and €2 = —1, this function reads, taking into account the representation ( |2.29| ) 
and introducing a dimensional regularization c? = 2 + 2 e because of the divergency when x^ — > 0, 



(A.51) 
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tUi Jo 
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r(l + e) Jo '^""'^ i<r''-''''''^''K^2^.. [q<^a,{l - -h (1 - 2a.) J , (A.52) 



28 



where we have performed the integration with respect to the angle {q,a/j). One can now perform 
the integration with respect to x'^, using the formula (6.576) of Ref. p^ : 

/ x-^KJax)J^ibx) dx=-i-] -] ' / \ ' 

/(T-A + l + pa-A + l- p }P'\ 



Here a = 0, p = —2iiyi, b = |(1 — 2ai), a = q^ai{l — ai) and —A = 2ieiP — 1 + 2e. Thus, 
On.i/. i + 2z/,; 1 1 + e 4 Jo V2/ 



l(l-2ai)2^ 



X TiieiU - iui + e)T(ieii^ + ii>i + e) 2-F1 icji/ - zi/j + e, ze^z/ + zz/j + e; 1; -- — . (A.54) 

\ Aai[l-ai)) 

One can now perform the integration with respect to z/j. We thus define 

9'{'y,x,,q) = j ^ ^El^'^^*{x,)UX^,q) exp (^^^;^(0, j . (A.55) 

The analytic structure of the integrand is very simple: it has poles at Vi = ±(ej u — ie). The 
integration contour can be closed either in the upper plane or in the lower plane, and give the 
same contribution. For we close the contour in the lower plane so as to pick up the residue 
at z/i = z/, and for g"^ we close the contour in the lower plane so as to pick up the residue at 
1^2 = —i^- Since 

/ l(l-2ai)2\ / l(l-2a2)'\ ^ 

2F1 0, 2tiy; 1; ^ ] = 2F1 [ -2zu, 0; 1; ^ = 1, A.56 

V 4ai(l-ai)y \ Aa2{l-a2)J 

one gets 

%,e,v r^(^ + ItiV) Jo ^ \ TT J 



(A.57) 



since x(0, u) is an even function of z/. Performing the integration with respect to 



1 1 . Y^(-ieu + -) 

rfa,(a,(l-a,))-2-- = -A-^^, (A.58) 
T[—2ieiV + 1) 



now reads 



Tial X1X2 f d'^q -iq.b r(2zz/) r(-2iz/) 

2r l^^2^iyj (2^'' ""V^o,-. r(l + 22z/) r(l-22z/) 

X i^J'^lxO i^J-'^lx^) exp (^^X(0, i^)Y^ , (A.59) 
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when taking into account that Y = Yi + Y2. Note that formula (|A.59|) could be equivalently 
obtained from Eq. ( [A.54| ) by performing the change of variable ai ^ Zi = {1 — 2aj)^ and 



using the Watson theorem |^3|. Taking the conjugate of formula (|2.29|) and then performing 
the change of variable a = a' — 1, one can verify that 



Using formula (|2.38|) , one finally gets 
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that is 



/ 7S^"'2V*(-.)Sr(.,)exp(^^,(0,.)F 
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r 2 , TT^ xlT^v{xi,X2,h,Y) = - \ xln^{x2,xi,b,Y). 



(A.60) 
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(A.62) 



4/ ^v- ' 4; 

which differs from the exact result ( 2.59| ) by a factor |. Thus, the approximation of using the 
elementary dipole-dipole total cross-section only gives half the correct result. 

One expects that the previous approximations should be correct when calculating the 



total cross-section. Indeed, integrating ( |A.45| ) with respect to b, one gets 



F^'\xi,X2,Y) 

/ X-^ X2 

From formulae ( p.65| ) and ( [A.14| ), this yields 
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(A.63) 
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Integrating with respect to x'^ and and using relation (p.67|) , the only remaining terms are 
rii = n2 = 0, due to the conservation of conformal weight, and one gets 

\ TT 



F^'\x,,X2,Y) 



2n 1^2 



I \l-2iv I \l+2iu 
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that is 
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{.X{.nu T^i)Yi + x('^2, T^2)Y2) 

(A.65) 

2Ho,v}{^ii^2), (A.66) 



which is exactly the result ( p.71| ). 



A. 6 Expression of E^'^ in terms of Bessel functions 



In this appendix we prove the expression ( p.30| ) of E^'^ . Starting from equation (|2.29|) and 
making the change of variable a = sin^ |, the mixed function E^'^ reads 



2iu 



TT Z=COSI 

e 2 K. 







2iu I ysint) dt. (A.67) 
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Let us calculate 
C 



e 2 K. 



2iv 



qp 



sin t] dt + 



.Q-P , 

e 2 K. 



2iv I Y sint 1 dt. 



Making the change of variable t' = vr — t in the second integral, one gets 



q.p 



C = 2 cos I =^ cos t K_2iu — sin t dt 



QP 



Combining 



and 



TT 



2 sin(7r/i) 



[/_,(A)-/,(A)] 



I,{\)=e-'^'^J,{e'^X), 
the Bessel function K can be written as 



(A.68) 



(A.69) 

(A.70) 
(A.71) 



K^iX) = -r(^)r(l - /i) [e'i^ J_^(e*5A) - e~^5^ J^(e^iA) 
C can then be computed using relation (6.688) of Ref. [p^] ) 
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(A.73) 



In our case /i = — 2ii^, a; = =2= and 2; = e*2 ^. Defining \E' the angle (g, p), and computing 



(A. 74) 



one gets 



c = -r(-2ii/) r(i + 2ii/) 



(A.75) 



where we have used the relation 



J, (e'^z) = e'^^J,{z) (A.76) 

in order to get rid of the minus sign arising from the ?/_ contribution. Combining expressions 
( |A.67| ) and ( [A.75| ), and using the expression ( p.22| ) for 6„ ,^, one finally gets 
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(A.77) 



in agreement with the more general result (10) of Ref. P2 
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A. 7 Properties of the three-points correlation functions E'^'" and 
q 

In this appendix we derive various useful mathematical formulae for the functions E'^''^ and 
E""'*^. Let us first show that E^''^ and E~'^~'^ are related by the following expression 



i^"''^*(pi0,P20) 
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which corrects formula (A. 12) of E^. Consider 
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(A.79) 



Using conformal invariance, we can take p2 — ^ C)0 in order to simplify the calculation and restore 
the p2 dependence afterwards by requiring the correct conformal transformation property. T 
then reads 



T 



1 



.Pic 



iPio'l IPoc 



'2+Aiu I Po'O 
\P00' 



1 



PlO' 



l+2ii/ 



\R' ~R\ 



-2+Aiv 



R'* - R* 
R'-R 



R' + l 



1+211/ 



(A.80) 



where pio = i? + |, p2o = -R — f , Pic = -R' + f and p20' = -R' ~ f • After performing the change 
of variable R = R' — R and introducing R = (r+^z,T reads 



where 
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(A.81) 



(A.82) 



Using the techniques developped in Ref. [Q, this integral can be computed after performing 
a Wick rotation for the y integration. The corresponding replacement y zye"^*^ reads 
z + z* = a + (3 and z — z* = {a — /?)e~^*^ with d?z = dxdy = ^dz dz* = ^da dj3. Separating the 
integration in a in three domains, the non-zero contribution is obtained for a G [—1, 0] due to 
the ie. Closing the integration contour for the (3 integration around the singularity /3 = 0, this 
yields 
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rf/?-^^-HTT— ,(A.83) 



where we have performed the change of variable a —a. These integrals lead to (3 functions. 
After some straightforward calculations, one gets 
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(A.84) 
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Combining formulae (|A . 78| , [Al8ll , |A .841) and using 



i^"''^*(pi0,P20) 



'R* + 



p \ 2 



(A.85) 



one finally gets formula (|A.78|) after restoring the correct dependence in p2- 

Note that from Eq. ( [A. 781 ) one can obtain the corresponding relation between E^''^* and 
E^^'*^. Indeed, performing the Fourier transform of both sides, one gets 
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(A.86) 



where the last integral is computed performing the integration (Pp^, = (^Pqiq with radial coor- 
dinates po'o = Ipo'oI exp(i0) (see the result (A. 11) of Ref. ||29[). Using the definition ( p.21|) . 



bn,u J \p\ V 2 2/ bn,u J \p\ V 2 2 

(A.87) 



where we have used the fact that 

i5"'''(pi,P2) = (-l)"i?"'^(p2,Pl) =i5;"'^(-p2,-pi). 

Thus, using formula ( |2.38D , one finally gets from Eq. (|A.86|) 



E^ip) = 1^1 



(A.88) 



(A.89) 



in agreement with formula (A. 15) of Ref. ||29|| . 



Note that the function E'"-''-' also possesses the following property 



Er*iP) 



27r2 r d^R -tq.R 



b* 



-e 'l-^E' 



\P\ 



R+P,R-P 
2' 2 



-R+P,-R-P 
2' 2 

where we have used the relation 

E-'-'*{p,,p,)=E-'-%-p,,-p,). 



E-^'-'ip), (A.90) 



(A.91) 
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Let us now prove the following orthonormalisation property for the three-points correla- 



tion functions E'^''^ 



j2 ' / \ 



(A.92) 



Using the definition ( ^.211 ), the left-handside reads 



X 



d^Rd^BL iq.{R-BL) pn. 



(A.93) 



Performing the changes of variables 



37^ 

Pi - po = ^ + y , P2 - Po = - y , Pi - Po' = i?' + y, P2 - Po' = - y , (A.94) 



and using the equality 

d^Rd^g d^x^ = d^p^ d^p^ d^p^„ 
one can apply the orthonormalization condition ( |2.48| ), which yields 



(A.95) 



d'x!, 

27TX[ 



27r2 27r2 1 



2^q 



h* h 2tt 



POO' 

vPoo' 



ni 



(A.96) 



Using the result (A. 11) of Ref. for the last integral, one finally gets the expected result 
(^). 



A. 8 Calculation of F^'^^ in the mixed representation 

In this appendix we show how to get the onium-onium scattering amplitude at fixed impact 
parameter using the mixed representation for the dipole distribution (see Eq. (|2.23|) ) and for 
the elementary dipole-dipole cross-section (see Eq. ( |A.3CI| )). Using the Fourier representation 
of the S distribution, formula (12.71) reads 



F«(xi,X2,F,&) 



r d^d^,2, ,2^ ,2(,. _ f -h- +g-b) 

2nx[^2TTx'2^ 



i2n) 



xn{xi, x'l , &i , ^1 ) n(x2, X2 , 62 5 ^2) o"dd (^'1 , ^2 , - ^2) ■ 
Combining this expression with formulae ( |2.23| ) and ( |A.30| ), one gets 

F«(xi,X2,F,6) - -i f -^^-'^-^ f ^'-^ ^ 



(A.97) 



J {2^Y ~ J 2nx'^2r^x',' ^^_J-oo 2^ x',^ ^"^^^""^ ^""^^ 



+00 

712 = — 00 ■ 



X2 E,-^'^^(x^) 2a^ "2 



27r 



x\x'^ ±2? /•+00 



16 



E 



4)' 



n— 1 



n+1 
2 



X EY\x\)EY{A) exp 



2a2iv, 



TT 



(A.9J 
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From the property (|A.9CI|) , one has 
condition ( [A.92| ), which yields 



n2,u2 _ n2, 1/2* (3g^j^ j-^Q^ apply the orthonormalization 



F^'\xi,X2,Y,b) 



iTa^^XiX2 f d?q -iq.b 



Y: / duE-/\x,)E-/{x,) 



2 I n-l\ ,.2 I /n+i 



exp 



TT 



{x{n,,u,)Y^ + x{n,y)%)]. (A.99) 



From the expression of the dipole density (|2.23|) , this finally reads 



i^{„;^}(a;i,X2,6) 



2 I / "-l 



^{n,i.}(^l,^2,^)- (A. 100) 



which is identical to Eq. ( p.58|) as expected. 
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